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Abstract 

An analysis is given as to why one can not directly evaluate continuum moment equations, i.e., 
equations involving powers of the position variable times charge, current, or energy/momentum 
operators, on the lattice. I examine two cases: a three point function evaluation of the nucleon 
magnetic moment and a four point function (charge overlap) evaluation of the pseudoscalar charge 
radius. 
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I. INTRODUCTION 



Beautiful hadronic moment equations can be derived from continuum field theory three 
and four point matrix elements |TJ. It is very tempting to try to implement these equations 
directly in lattice simulations. A number of papers in fact have used them, assumed them to 
be true, or examined their consequences in lattice calculations^, ||, f|, ||. Physical quantities 
considered have been charge radii[^|, ||, magnetic moments^, §] and quark total angular 
momentum [[J . However, these type of equations share one crucial feature in their derivation: 
a derivative with respect to momentum transfer evaluated at zero momentum. This last step 
can not be reproduced on the lattice because of the finite momentums available there, so 
the question arises as to the validity of such continuum-derived expressions evalauted on the 
lattice. We will examine two such expressions in this paper and will see that the continuum 
expectations and the lattice reality can differ markedly. 

II. TWO EXAMPLES 

Let us recap the situation for one such specific case, an expression for nucleon magnetic 
moments given in Refs.p, ||. Both two and three point functions appear in this expres- 
sion. The time ordered two point function, using the proton interpolation field, Xa( x )i is 
(understood a, a' sums) 

G pp (t;p,T) = Y.^ ip ^ a > a (™ C \T{xi{x)x P ai m |vac>, (1) 

X 

t 




and Tfc = i 






expression for r = I^ is given by 
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The long Euclidean time limit of this 



M *' P,l4j 2E p (2«)3 6 ' {I) 

where is the (dimensionless) nucleon mass, "a" is the lattice spacing, and k = 
l/(2(m + 4)), In addition, Z is the normalization factor, (vac|x„ ;cori< (0)|p, s) = Zu a (p,s), 
where u a (p, s) is the free Dirac spinor, \p, s) is a continuum proton state, and Xa Cont ( x ) is the 
continuum interpolation field. The three point function we need, which uses the conserved 
vector current, Jfi(x), is 

G pJ ^(t 2 ,t i; p,f,T)^-tY,e-^ 2 e^ (3) 

X2,Xl 



which has the long Euclidean time limit, 

G pJjP (h,t i; 0,-q,T k ) fo-*^ 1 _LJ^ e -» W (fe-<i) e -^*iey fc (^) 4 G ra (g2). (4) 

The Minkowski four momentum transfer squared is given by Q\ = 2m^{Ei — mjv), we are 
assuming continuum dispersion, E\ = m 2 N + qf , and (<^)j = ^, £ = 0, ±1, ±2, . . . , N in 
a given momentum direction for a square spatial lattice of size N s = (2N) 3 . (All these 
considerations can be reformulated for a lattice with an odd number of spatial sites in 
given directions, but no fundamentally different conclusions or observations results.) Taking 
a "continuum" derivative of G p j. p (t 2 , t\\ 0, — q, T k ) with respect to (qe)i, evaluated at zero 
momentum, gives us another three point function, which we will define as 

G pJjP {t 2 M] (xi)i,T k ) = ( r fe)a'a( vac IXa(a;2)(^i)i^(^i)Xa'(0)|vac). (5) 

X2,Xl 

Setting this quantity equal to the continuum derivative of Eq. (H) then results in 



G p j jP (t2, ti] (xi)i, T k ) 



G pp (t2',p, T4) 



(*2-*i),ti>i G m (0) 



D 



The "D" notation reminds us that this result follows from using a momentum derivative 
evaluated at zero momentum. 

Ref.|Q justified this procedure in the following way. One can imagine first taking the 
spatial momentum derivative of the continuum analog of G p j, p (t2, t\\ 0, — q, T/.), evaluated 
at q = 0, and dividing by the continuum two point function. One then transcribes this 
result into lattice language by changing the continuum matrix elements into lattice ones and 
making the appropriate substitutions for the spatial integration and the various fields. It 
was found there that this procedure resulted in a lattice measurement giving unrealistically 
small neutron and proton magnetic moments, G m (0). In particular, there was a downward 
trend in the data for smaller quark mass. The results in are similar, although the values 
are larger. 

In order to understand why this equation fails on the lattice, let us expand the position 
variable in terms of the momentum eigenstates of the periodic lattice. In one dimension we 
have the discrete completeness statement 

1 N 

— E e-*» = <&,„. (7) 

ZJV „=_jV+l 
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On a periodic lattice with an even number of sites, choosing an origin forces one side or 
the other of the lattice to have one extra spatial site. For this purpose, one needs to define 
a function which represents a linear function everywhere except at the extra site, n = N, 
where, because it may be considered equally distant from the origin in either direction, we 
will take it to be zero. (This is the same function used in the numerical evaluations in 
Thus, 

n , n 7^ N 
,n = N 

We expand this in terms of the momentum eigenf unctions, e mn , 



F{n) 



(8) 



N 

F{n)= Yl Cze iqen . 

l=~N+l 



Using Eq.(^), this gives 



1 



N-l 



C t = y F(n) e 

2N ^ y ' 



■iqen 



(9) 



(10) 



n=-N+l 

Summing this finite series in the usual way by multiplying both sides by a phase factor, 
e~ m , and shifting the summation limits (note that e ±mN = (—1)^), one finds that 



C, 



This leads to 



Gpp(t2',P, T4) 



(t 2 -ti),ti >i 



i(-l)^cot(©/2) ,£^0 
J = 



p m N zi p 

-e ijk J2(-lY cot(q e /2) - 



(11) 



Am 



s ^"" N 1^0 

The "S" notation now reminds us that this result follows from explicitly performing the lat- 
tice sum. The coordinate when inserted in Eq.(|5]), projects over the lattice momentums 
with a function given by Eq. ( [11]) . 

The leading terms in Eq.flT^D define what we will term the extreme Euclidean time limit 
(EETL) on the t\ variable as, 

G p j jP (t2, t\] (xi)i, Tfc) 

s 



E f 



-{qt)iG m {Ql). (12) 



„(m,]v— -Ei) ti 

eetl e ____r (n^ 



(13) 



Gpp(t2', p, r 4 ) 

where we have approximated sin((<fi)j/2) ~ (ft.)i/2. Eq.([T3|) is not suitable to measure the 
magnetic moment on the lattice. At fixed finite (<fi)i, the signal involves only the lowest 
nonzero component of G m (Qf), and is not time independent. 
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FIG. 1: The magnetic moment of the proton at four values of m q , extracted three different ways 
on a 16 3 x 24 lattice from the data of Ref.[§]. Square symbols are extrapolated zero momentum 
form factors. The diamonds are extracted values assuming Eq.(^) were valid. Circles show the 
result one would find for Eq.(^) if the lowest momentum EETL limit in Eq.(|i~3|) was dominant in 
the momentum sum in Eq.([l2]). 

It is not a contradiction that Eqs.(|^) and ([13]) disagree with one another even in the 
qi — > (or N — > oo) limit. If one formed a discrete lattice derivative for the left side of Eq.(^), 
equivalent to simply evaluating Eq.([|) at the lowest spatial momentum and dividing by that 
momentum, one would obtain a result consistent with (|6|) in the q\ — > limit. Of course, 
one may always use external field methods to consistently extract magnetic moments. Then 
one is effectively taking derivatives with respect to the external field rather than momentum 
to isolate the coupling. 

Fig. 1 shows measurements of the lattice proton magnetic form factor at zero momentum 
transfer on a 16 3 x 24 lattice. The open square symbols are extrapolated from nonzero 
momemtumflJ. A measurement assuming Eq.(§) yields the solid diamond symbols, which 
are trending downward as a function of decreasing quark mass. The results in Ref.Q are 
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similar. Using Eq. (|T3| ) one can test the extent to which the lowest momentum is contributing 
to the result from fll^D because one has separate data for G m (Qf) at the lowest momentum 
transfer^]. The solid circles show the result one would find for Eq.(||) if the lowest mo- 
mentum EETL limit in Eq.(|T^) was dominant in the momentum sum. The trend is also 
downward for smaller quark mass, which is a result of the decreasing value of G m (Ql) as well 
as the greater exponential suppression from the e*™^ - ^ 4 factor, (t = 7~ for the data from 
Ref.|Q). The qualitative behaviors are remarkably similar for these two Eq.@ measure- 
ments, although significant cancellation is probably occurring for the diamond data from 
the {-l) e factor in Eq.©. It was speculated in [[3], f|] that the downward trend possibly 
reflected the fact that the nucleon was not well contained in the lattice volume. The Fig. 
1 data strongly suggests that the downward trend in the diamond symbols is instead the 
result of two factors, decreasing G m (Qj) and greater e^ riN ~ E ' t suppression, similar to the 
EETL case[|. 

As another example of a continuum moment equation on the lattice, consider the charge 
overlap measurement of the pion charge radius. Time separated measurements of charge 
overlap matrix elements were first considered in lattice calculations in Ref.||. 

We will start with the results derived in Ref[0. Using u, d flavor conserved lattice charge 
densities, p u,d (x), one has for t 3 ^> ti j2 3> 1, 

Ea, * < 0\tf(x 3 )T (p u (x 2 )p d (x 1 )) <f>{z)\0 > 



E^<o|0t(x 3 )<Kz)|o> >v: (14) 

where z = (z,0) and is a charged pion interpolation field. V% d (r,t) can be written 
using space and time translational invariance between the zero momentum pion states (r = 
x 2 — xi, t = t 2 — ti) as 

V: d (r,t) =< n(0)\T(p u (r)p d (0))\n(0) > . (15) 
The discrete Fourier transform is defined to be, 

Q(g 2 , t)=N s Y; e-^ r V: d (r, t). (16) 

r 

It is important to let t 3> 1 in Euclidean space in order to damp out the contributions of 
higher mass intermediate states when a complete set of states is inserted between the charge 
densities in (|l5l). One then has 

Q(q\ t) ^ (^ + ^)^ (Q 2 )e^-^ (17) 



where Q 2 = 2m n (E q — m n ) and F n (Q 2 ) is the pion form factor. 

A continuum derivative of Eq. (|T6| ) with respect to q 2 at zero momentum forms the quan- 
tity 

TZl(t) = N s J2r 2 V: d (f,t). (18) 

r 

Using Eqs. (p!6|) , (P~^|) and ([H]) and following a procedure similar to the above for the magnetic 
moment, we obtain the continuum derivative result, 

nl(t)\ D ^2Rl d + ^, (19) 

where R 2 d is the charged pion u, d quark charge radius. This is the same result as in Ref. || 
when the bag sources there are replaced with zero momentum pion sources. One could 
imagine evaluating this expression on the lattice to try to extract the charge radius from 
the time constant term, but we will see this hope is ill-founded. 

In order to explain why one can not measure the charge radius from an expression like 
Eq. (|I8"D , let us now expand the square of the position variable, r 2 , in terms of the momentum 
eigenstates. For a one dimensional lattice we will consider, 



N 

9 

n 

e=-N+i 



2 - Yl K e e iqin . (20) 



We then find that the coefficients, Kt, are given by 

|(-l)Wfe/2) ,^0 

'i(iv 2 + i) ,e = o 

In a three dimensional context, the quantities on the right in ( |2"T| ) are multiplied by zero 
momentum Kronecker deltas in the transverse directions. 

Using Eqs.(|2H|) and ( PT|) in ([TB| ) and inserting a complete set, one finds that 

TZl(t)\ m (N 2 + \) + I £(-1)' csc 2 (W2) {E " + m / F 2 M^- E ^. (22) 



We can define a discrete lattice charge radius as 



In the EETL, we then find, 



- 3 (1 F l m \ (23) 



K(t)\ s ™ (N 2 + 1)- 12(^ - \{Rl 4 ) L )e^-^. (24) 



So we see that in contrast to Eq.(|T9D, the first (time independent) term is essentially mean- 
ingless. It is conceivable that one could extract a measurement of (R^ d ) l from the second, 
time dependent term, but it would be much simpler to project out the form factor at the 
lowest lattice momentum in the usual manner and use Eq.(^3]). Again, the fact that ( |T9"D 



and ( pip do not agree even for very fine lattices (N — > oo) is not a contradiction because 
one has not taken a discrete lattice derivative, but a continuum one in producing (|Hf). 

III. SUMMARY 

We have seen why it is not possible to directly evaluate continuum moment equations 
on a periodic lattice. Continuum moments of lattice operators in a periodic system do not 
project onto good momentum and so do not isolate low momentum properties. The present 
author pointed this fact out some years ago in the context of charged pion polarizability 
calculations 0. In a more general sense, the lesson we have learned is that in order to 
deduce continuum properties from the lattice, it is important to treat the lattice as a self- 
consistent physical system. Position functions are meaningful only if expanded in terms of 
the available momentum eigenstates of the system. It can be very misleading in general 
to try to take continuum field theory equations and simply "latticize" them in deducing 
physical properties. 
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